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We derive two-loop renormalization-group equations for tiie lialf-filled one-dimensional Hubbard chains cou- 
pled by the interchain hopping. Our renormalization-group scheme for the quasi-one-dimensional electron sys- 
tem is a natural extension of that for the purely one-dimensional systems in the sense that transverse-momentum 
dependences are introduced in the g-ological coupling constants and we regard the transverse momentum as a 
patch index. We develop symmetry arguments for the particle-hole symmetric half-filled Hubbard model and 
obtain constraints on the g-ological coupling constants by which resultant renormalization equations are given 
in a compact form. By solving the renormalization-group equations numerically, we estimate the magnitude of 
excitation gaps and clarify that the charge gap is suppressed due to the interchain hopping but is always finite 
even for the relevant interchain hopping. To show the validity of the present analysis, we also apply this to 
the two-leg ladder system. By utilizing the field-theoretical bosonization and fermionization method, we derive 
low-energy effective theory and analyze the magnitude of all the excitation gaps in detail. It is shown that the 
low-energy excitations in the two-leg Hubbard ladder have SO(3)xSO(3)xU(l) symmetry when the interchain 
hopping exceeds the magnitude of the charge gap. 
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I. INTRODUCTION 

Renormalization-group (RG) method is one of the most 
powerful and promising tools to tackle low-dimensional elec- 
tron and spin systems.' It has a long history of research es- 
pecially on one dimensional (ID) systems, since the RG the- 
ory is superior to take into account low-dimensional compet- 
ing fluctuation effects, i.e., it can sum up systematically the 
logarithmic-singular particle-particle and particle-hole chan- 
nels which appear in all order of perturbation theory. Siii 
It has been clarified that the RG method describes various 
ID ground states: the Tomonaga-Luttinger (TL) liquid state, 
the charge-gapped Mott insulating state at half filling, and 
also the spin-gapped Luther-Emery state. ''^ Not only for the 
most divergent terms, the next-to-leading logarithmic singu- 
lar terms have also been studied based on the two-loop for- 
mulation of the RG theory,^ ** where singular self-energy cor- 
rections in addition to the vertex corrections are taken into 
account. Recently the RG theory is generalized to apply to 
two-dimensional electron systems.'' The main difficulty in the 
RG formulation for two-dimensional systems resides in the 
fact that the momentum dependence of the coupling con- 
stants is essential but the number of independent coupling 
constants is large and it becomes hard to analyze the RG equa- 
tions even for the one-loop level. Several attempts have been 
made by focusing only on dominant scattering processes on 
the Fermi surface^ and by discretizing the Fermi surface into 
finite number of pieces, i.e., so-called patches.** For electron 
systems in arbitrary dimension, a nonperturbative RG theory 
has also been formulated"^ and has been applied to two dimen- 
sional electron systems by considering leading two particle 
interactions, i.e., within the one-loop level. " Quite recently 
the effect of the two-loop self-energy corrections have been 
examined, 'iiSiiiiM. while the two-loop vertex corrections are 
considered only for the system with flat Fermi surfacepi^iii 

In quasi-one-dimensional (QID) electron systems, the im- 
portant issue to be clarified is the dimensional crossover from 



one to higher dimensions which would occur by changing pa- 
rameters or temperature .-i-^"* In real QID compounds, the TL- 
liquid behavior is expected at high temperature, however, the 
effect of warping of the Fermi surface due to the small but fi- 
nite interchain hopping is enhanced at low temperature where 
the Fermi-liquid behavior can be expected if the system is 
metallic, and finally the system has an instability to symmetry- 
broken states. The RG approach is also powerful and succeeds 
in the description of these physical picturespifii^iiS In the early 
RG analysis, the effect of the one-particle interchain hopping 
is treated perturbatively, however, it is found to be relevant 
even in the noninteracting case and the perturbative treatment 
is invalid at low temperature. In order to clarify the dimen- 
sional crossover phenomena properly, one has to formulate the 
RG with the nonperturbative treatment of the interchain hop- 
ping, i.e., based on the warped Fermi surface. In this sense, the 
formulation is analogous to that in the two-dimensional RG 
scheme since one has to discretize the Fermi surface. In the 
QID case, the RG has been formulated by considering finite 
number of chains N± (Af_L-chain RG scheme)22i^ii22i2i2iS^ 
where the transverse momentum is regarded as a patch in- 
dex. Based on this scheme, the QID systems have been an- 
alyzed intensively within the one-loop ievel22i^i*2SiSi21i2^ and 
the self-energy corrections have also been investigated.^^ At 
commensurate band filling, the dimensional crossover prob- 
lem becomes nontrivial since the electronic correlation has 
the strongest effect and leads to the Mott insulating state if 
the system is half filled. The effect of the umklapp scattering 
between electrons, which is a trigger for the ID Mott insu- 
lator, has been investigated by the one-loop RG,^^ however, 
in order to clarify the electronic states in the Mott insulator 
one has to examine the properties of the one-particle Green's 
function, i.e., the self-energy corrections, whose singular con- 
tributions only appear beyond the one-loop level. The effects 
of the two-loop self-energy corrections have also been exam- 
ined in the QID systems without considering two-loop vertex 
corrections^^ however, a systematic two-loop RG including 
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both the two-loop vertex and self-energy corrections is not 
formulated yet. Recently this issue has also been addressed 
by a numerical method expanding the dynamical mean-field 
approach (chain-DMFT)2i22i22i^ and by a field-theoretical 
method with the RPA treatment of the interchain hopping.^' 

From a technical point of view, it is generally hard to 
gain physical insights of results of scaling flows in the QID 
RG, since the number of independent coupling constants be- 
comes large as N± increases. As a minimal system of the 
coupled chains, one can consider a two-leg ladder system 
{N± — 2). The two-leg ladder system itself has nontriv- 
ial and interesting features^- and has been examined inten- 
sively by using the RG method^^'?^i^^'?^i"i?^'?Wi^^ and also 
by the high-accuracy numerical technique called the density- 
matrix-renormalization-group (DMRG) method,'*^ '*'* where it 
has been confirmed that both the charge and spin modes have 
excitation gaps for the half-filled Hubbard ladder. In this anal- 
ysis, one can easily see that a naive one-loop RG analysis of 
the excitation gaps is not satisfactory since the RG method 
breaks down at a energy scale corresponding to the largest ex- 
citation gap in a system. In order to analyze the lower-energy 
properties, one has to derive an effective theory by tracing 
out the gapped modes based on the field-theoretical bosoniza- 
tion/fermionization treatment. As for the two-leg Hubbard 
ladder, Lin, Balents and Fisher^'^ obtained the SO(8) Gross- 
Neveu model as an effective theory in the low-energy limit 
and examined the excitation spectrum. The extended two- 
leg Hubbard model including additional interactions is also 
examined"^" '*' and quantum phase transitions between com- 
peting ground states have been clarified in this context of the 
one-loop RG. Despite that the analysis of the two-leg lad- 
der systems based on the one-loop RG succeeds in describ- 
ing the ground state properties, it is not easy to extend the 
analysis to the case with large number of chains, since the 
field-theoretical approach is restricted to the small number of 
chains. In order to overcome this problem, we formulate, in 
the present paper, the two-loop RG theory for the QID elec- 
tron systems. Even in the two-loop level, the perturbative 
approach also breaks down at energy scales of the excitation 
gaps, however, the respective excitation gaps can be estimated 
by analyzing the scaling behavior of the couplings for respec- 
tive modes, without following the tracing-out procedure. We 
confirm that the present scheme works even if the respective 
modes are not independent by revisiting the two-leg ladder 
systems. 

This paper is organized as follows. In Sec. HI] we introduce 
the finite A^x -chain half-filled Hubbard model coupled by the 
one-particle interchain hopping, and derive the corresponding 
5-ology model by linearizing the energy dispersion where the 
effect of the interchain hopping is treated nonperturbatively. 
By developing symmetry arguments for the particle-hole sym- 
metric half-filled Hubbard model, we obtain constraints on the 
g-ological coupling constants. In Sec. |lll] we formulate the 
RG based on the Kadanoff-Wilson approach up to the two- 
loop level, where vertex corrections are taken into account 
based on the third-order perturbation theory, in addition to the 
second-order calculation for the self-energy corrections. Re- 
flecting the symmetries that the particle-hole symmetric Hub- 



bard model has, the resultant RG equations can be written in 
a compact form where the physical picture can easily be cap- 
tured. By solving the RG equations numerically, we estimate 
the magnitude of the charge and spin gaps. In Sec. lIVI in or- 
der to indicate the validity of the present method, we consider 
a most simple but nontrivial case N± = 2, which corresponds 
to the two-leg ladder, and analyze the excitation properties 
in detail by combining the field-theoretical bosonization and 
fermionization method. Finally, the results are summarized in 
Sec. |V] Technical details are given in the Appendices A and 
B. In the Appendix C, we give a related issue which supports 
strongly the validity of the present estimation of excitation 
gaps. 



II. MODEL AND SYMMETRY ARGUMENTS 

We consider the bipartite QID Hubbard model at half filling 
with t|| ^ t±, where the transfer integral along chains is i|| 
and that between chains is t±. Our Hamiltonian is given by 

+ C/^nj,;,tnj,;,x, (2.1) 

3,1 

where ; ^ is the annihilation operator of electron on the jth 
site in the l\h chain with spin s, and Uj^i^s = j ^c^ i s ~ \- 
The system size along chains ( A^n ) is considered to be suffi- 
ciently large and the sum of the site index, which runs j = 
1, • • • , Af||, is to be understood as an integral in the thermody- 
namic limit. The chain index runs / = 1 , • • • , iVj^ and we con- 
sider the system with finite number of chains iVj^ where the 
periodic boundary condition is imposed Cj,N^+i,s = Cj,i,s- 




FIG. 1: (Color online) Fermi points (closed circles) in the present 
half-filled QID Hubbard model with the periodic boundary condition 
in the transverse direction. The case for N± = 8 is shown. 
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A. gr-ology notation 
The kinetic term of the Hamiltonian is given by 
Ho = ^£(fe)ct(fc)c,(fc), 

e{k) = ^11 ^ 



(2.2) 



(2.3) 



where fe — (fc||, fc±) and the lattice constant is set to unity. 
Since the system is particle-hole symmetric, we can assume 
t|| > and t± > without losing generality. Since the num- 
ber of chains N± is finite, the transverse momentum is given 
by 



k I = n. 
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(2.4) 



where [x] is the Gauss symbol denoting maximum integer 
which does not exceed x. By assuming t± ^ t^^, we lin- 
earize the dispersion where the situation can be simplified as 
follows. Up to the lowest order in t± the kinetic term with the 
linearized dispersion is given by 

Ho = ^ep(fc)ct^,(fc)cp,,(fe), (2.5) 

k,p,s 

£p{k) = v{pk\\ — kp) — 2t± COS k±, (2.6) 

where v = 2t|| and hp = 7r/2. We introduce the bandwidth 
cutoff A. In this approximation, the warped open Fermi sur- 
face (Fig.[0 is specified as a function of k±: 



fc_F(fc±) = kp + 2— cos k±, 

V 



(2.7) 



and the energy dispersion J2.6t can be reexpressed as ep(fe) = 
— kp{k±)]. Thus we regard the transverse momen- 
tum k± as a patch index in the present RG formulation. The 
greatest merit of the present formulation lies in the fact that 
the transverse momentum k± is a conserved quantity, i.e., the 
patch index is a good quantum number and the ambiguity of 
selecting patch index disappears. 



Following the conventional y-ology approach,"* we clas- 
sify the interaction part of the Hamiltonian H\ = 
U J2j I into the forward, backward, and umklapp 

scattering processes, by focusing on the longitudinal momen- 
tum We introduce the coupling constants gi±, g2±, g\\, 
g3±, and g^, which represent the backward scattering with 
the opposite spins (gi±), the forward scattering with the op- 
posite spins {g2±), the forward scattering with the same spins 
(.g||), the umklapp scattering with the opposite spins ig3±), and 
the umklapp scattering with the same spins (gay). In terms of 
the Hubbard interaction U, the magnitudes of the couplings 
are given by gi± = g2± g3± = U and .g[| = 53 1| = 0. The 
.(7|| and g-^ processes are absent in the original Hubbard inter- 
actions, however, can become finite under the RG scaling pro- 
cedure. Furthermore, the coupling constants are differently 
renormalized depending on the external transverse momenta 
of the vertex and have the explicit transverse-momentum (i.e., 
patch-index) dependence. To take into account these effects, 
we formally introduce the transverse momentum dependence 
of the coupling constants in the initial g-ology Hamiltonian. 
In the most general form, the interaction part of the Hamilto- 
nian is given by 



Hi 



fci ,k2,q,s 



+ 92Uq±:k^i,k^.,) C4.,,(fc2) C^-,s{k2 - Q) C_^s{ki - Q) 



\ H S'||(qx,fcii,fci2) c+^,(fci) c+,^(A;2) c^-A^i - Q) c^Jki - Q) 



fcl ,k2,q,s 



2V 
1 

2V 



cljki) c_Jki - Q) cl^,{k2) c_-Ak2 +Q-G)+ H.c. 



fcl ,fc2 ,g,s 



Y 9mq±-k±i.k^2) c+^,(fci)c__,(fci-Q) ct^^,(A;2)c_^^(fc2 +H.C. 

ki,k2,q,s 



(2.8) 



where s =]{[) for s =i(t), and Q = (tt + q\\,q\_), G = 
{27T, 0), and V — N\\N±. The momenta fc||i and k\^2 are as- 
sumed to take values near kp{k±). In the transverse direction, 
on the other hand, the momenta k±i and q± can take the val- 
ues in — TT < k±i, q± < tt, and the momentum {k±i ± q±) is 



assumed to reduce the first Brillouin zone, then all the pos- 
sible scattering processes are taken into account, including 
the transverse umklapp scattering. The respective scattering 
processes are shown in Fig. |2] We will neglect the forward 
scattering with the same branch, so-called 54 term, since this 
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fei — Q, s 




ki — Q,s 
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+ Q-G,s 


k2,S 


fe2 + Q - G 



k2 — Q,S k2,S k2 — Q,s ki — Q,s k^ — Q^s ki — Q,s 

FIG. 2: g-ology notation. The solid (daslied) line denotes a right-moving (left-moving) electron, ki = (fc||i, fcxi), <? = (ti" + , fZi), and 
G = (27r,0). 



process does not show the logarithmic-singular behavior in 
perturbation and is known to yield only quantitative changes 
in velocities for the ID case. In terms of the Hubbard in- 
teraction U, the magnitudes of the couplings are given by 

9l±{q±,k^l,k^2) = 52±(gj_,fcii,fcj_2) = 93±{q±,k^i,k^2) = ^ 

and = 93Uq^,k^,,k^2} = 0- To simplify the 

notation, we will suppress the _L index of the transverse mo- 
mentum in the following. All the coupling constants are as- 
sumed to be real. In order to make Hi hermitian, the coupling 
constants must satisfy 



9l±{q,ki,k2) 
92±(q,ki,k2) 
9\\(qM,k2) 
93±{q,ki,k2) 
93\\{qM,k2) 



9l±{qM,ki), 

92±{q,k2,ki)> 

3||(g,fc2,fcl): 

93±{-q,k2,ki)i 

93\\(-q,k2,ki)- 



(2.9) 



As in the ID case, the physical picture becomes transparent 
by introducing a new set of the couplings: 



9p{qM,k2) = 92±{q.kiM) + 9\\{qM,k2)^ 

9a{qMM) = g2^{q,klM) ~ 9\\[q.M,k2)-' 

9c(qM,k2) = .93±(g,fci,7r-fe2)i (2-10) 

9s(q,kiM) = 5l±(g,fci,fc2)i 

9cs{q,ki,k2) = 53|| (g,fci ,7r — fc2) J 



where gp and gc (go- and gs) are the coupling constants repre- 
senting the charge (spin) degrees of freedom. This picture can 
easily be captured by noting that, if we neglect the momentum 
dependence of the coupling constants, the gi±, g2±, g\\, and 
g3± terms of the Hamiltonian ( 12. 8> are written in symmetric 
forms as'*^ 



where the respective chiral density operators are given by 

k.s 
k.s 



(2.12a) 
(2.12b) 



(2.12c) 



Jp'iQ) = $]cp,T(fe)cp,i((7r,^) -fc + q), (2.12d) 



and J+{q) = [Jp-(-q)]t, j;+iq) = In the ID 

half-filled Hubbard model (gp — g^ and 5^ = 9s), it is known 
that the charge part, in addition to the spin one, also becomes 
SU(2) symmetrici^ Even in the QID case, the model has an 
additional SU(2) symmetry, which is shown explicitly in Sec. 
Ill B 31 The Gcs coupling represents the spin-charge coupling 
term in the ID case as seen from the bosonization technique.— 
In the notation of Eq. ( I2.10> . the conditions of the hermitian 
(jZOj can be expressed as g^iqMM) = 9u(qM,ki)' for ^ = 
p, a, s and g^^q^k^M) = 9^(-q.,T^-k2,^^ki) for 1^ = 0, cs. The 
number of independent coupling constants gi in Eq. ( 12. 8> is 
[{N^ + l)/2. 



B. Symmetry arguments 

The Hubbard model ( 12. H is known to have high symme- 
tries, however, the g-ology Hamiltonian MM is generalized 
one including low symmetry. Reflecting symmetries that the 
Hubbard model has, there appear several constraints on the 
g-ological couplings and the resultant RG equations can be 
simplified. In this subsection, we clarify relations for the cou- 
pling constants protected by the symmetries. 



1. Spin-rotational SU (2) 

^ P'*^ The Hubbard model ( 12. Il l is invariant under spin-rotation, 

+ j;^(g)Ji^p(-q) + ^^ j;+(g)Ji-(-q),(2.11) while the .g-ology Hamiltonian ^ includes the spin- 

^ p q ^ p q anisotropic case. The spin-rotational symmetry can be argued 
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in terms of the generators of the spin rotation which are noth- 
ing but the spin operator: 



2 ^ 

fc.Si ,S2 



(2.13) 



The arbitrary global spin rotation by these generators can be 
represented by the SU(2) matrix: 



ha 

a* 



(2.14) 



where and ba are complex numbers satisfying la^P + 
I^CTp — 1- Obviously the Hubbard Hamiltonian i2. li is invari- 
ant under the transformation ( I2.14l i. By requiring the g-ology 
Hamiltonian ( 12. 8> to be invariant under this rotation, we ob- 
tain the constraints on the coupling constants. In the notation 
( I2.10> . the constraint relations are given by 

gs(qMM) ^ 9a(qMM)^ (2.15a) 
9c{q.ki,k2) ~ 9c{-!T — q+ki + k2.ki,k2) 

= 9cs(q.kiM) ~ 9cs(Tr-q+ki + k2.ki.k2)- (2.15b) 

Since the property of the spin-rotational invariance is hold un- 
der the RG procedure, these relations can be considered as 
constraints on the renormalized coupling constants. 



where the so-called 77-pairing operator is given by 
V = ? ^ct(fc)c4((7r,7r) - k). 



k 



(2.19) 



The arbitrary rotation by these generators can be represented 
by the SU(2) matrix: 



(2.20) 



where Op and bp are complex numbers satisfying \ap 



l,andzj./ = (— 1)-'+'. This transformation commutes 



with Eq. ( 12.141 1. One easily finds that this symmetry breaks 
down if the Hubbard model is extended, e.g., by including an 
additional intersite interaction. In the Fourier space with the 
linearized dispersion, the transformation ( I2.20t corresponds 
to 

Cp,-\(k) apCpj(fe) + 6pc|, j^((p7r, tt) - fc), (2.21a) 



ljk)~b;Cj,^^{{pTr,7r)-k). (2.21b) 



The kinetic term ( 12. 5t is invariant under this transformation. 
By imposing the condition that the g-ology Hamiltonian il.Ei 
is invariant under this transformation, we obtain 



2. Particle-hole symmetry 

The present bipartite half-filled system is invariant under 
the particle-hole transformation Cs{k) ^ ct((7r,7r) — fc), 
where Cg{k) is the Fourier transform of ; ^. In the linearized 
dispersion ( 12. 6> . this particle-hole transformation corresponds 
to 



Pi -5 



(2.16) 



In order to make this particle-hole symmetry meaningful, the 
number of chains N±_ must be even, otherwise the k±_ [Eq. 
MAW cannot become symmetric in this transformation. By 
imposing the condition that the y-ology Hamiltonian MM is 
invariant under this rotation, we obtain the constraints, in the 
notation ( I2.10L 

9y(qMiM) ~ 9y(-q,TT-ki,TT-k2)i (2.17) 

where i' — p,a, c, s, cs. We note that, by combining the rela- 
tion .gc/cs(g,fei,fc2) = 9c/cs{-q,-R~k2,7T-ki) [obtained from Eq. 
(|Z9}], we find 9c/cs{qMM) = 9c/cs{qMM)- 



3. Pseudospin SU(2) 

In addition to the particle-hole symmetry, the system has an 
additional symmetry, if the interaction is on-site one only."*^ 
The generators of this SU(2) are given by'*'' 



2i 



(2.18a) 

Q^^i^:ct(fc)c,(fc):, (2.18b) 



9c(q,ki,k2) ~l" 9c(iT — q+ki+k2,ki,k2) 

= +9p(qMM) + 9p(-K-q+ki+k2MM)^ (2.22a) 

9c(q,kx,k2) 9c(-n — q+ki+k2,ki,k2) 

= —9a(q,kiM) + 5o-(7r-g+fci+fc2,fcl,fc2)' (2.22b) 
9s{q,ki,k2) 9s{7r — q+ki+k2,ki,k2) 

= ~9cs{q.ki,k2) + 9cs(lT-q+ki+k2,kiM)- (2.22c) 

The first relation is a natural extension to the known relation 
for the purely ID case.'*^ The last two relations, which do not 
appear in the ID limit, imply that the couplings gc and {gs 
and gcs) are not independent and related to each other 

The relations M.22\ can also be derived from the spin SU(2) 
relations ( I2.15> by using the charge-spin duality relation, as 
explicitly shown in the AppendixIXl 



C. Two-loop RG theory for the ID Hubbard model 

We briefly recall the known results of the two-loop RG the- 
ory for the purely ID case, by focusing on the half-filled ID 
Hubbard model: 

i?iD = -^E +H-C.)+f/E"j'T"ja> (2-23) 



k,s 



where Cj s is the annihilation operator of electron on the jth 
site with spin s, and nj^g = Cj ^Cj ^ — ^- The linearized dis- 
persion is e{k) — — 2t cos k v{±k\^ — kp) where the Fermi 
velocity and the Fermi momentum are w = 2i and kp = 7r/2. 
The g-ological scattering matrices are the same as Fig. |2] and 
we introduce gp = {g2i_ + g\\), 9a = (52_l - 9c = 93±, 
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cjs = gi±, and Qcs = g3\\, as before. The two-loop RG equa- 
tions for the respective couplings are given hy^ 



jr" - 


1 O/^"^ 0/^ 


(2.24a) 


jr^ - 


+2GpGc — G^pGc — G^, 


(2.24b) 


jr^ - 


— 2Gg — '2GfjGg^ 


(2.24c) 


^G - 

jr- ~ 




(2.24d) 



where I is the scaling parameter and the initial values are given 
by Gi{0) — gi/{2TTv). We have neglected the Gcs coupling, 
since this has an irrelevant canonical dimension."*^ 

These RG equations can be simplified reflecting the sym- 
metries of the system. The spin-rotational SU(2) symmetry 
ensures Ga{l) = Gs{l), which is obtained from Eq. ( I2.15a> 
by neglecting the transverse momentum dependences. This 
relation holds even under the scaling procedure. The particle- 
hole symmetric Hubbard model has another pseudospin SU(2) 
symmetry and then the total Hamiltonian is characterized by 
the SU(2)xSU(2) symmetry. This pseudospin SU(2) sym- 
metry ensures Gp(/) — Gc{l), which can be obtained from 
Eq. ( l2.22aL and thus this can be considered as the "charge" 
SU(2) symmetry. In this SU(2)x SU(2) symmetric case, the 
RG equations J2.24t can be simplified as 

^Gp = +2G2-2G3, (2.25a) 

^G, = ~2Gl-2Gl, (2.25b) 
dl 

where the initial values are given by Gp(0) = Go-(O) = 
U /{2-Kv). For repulsive interaction J7 > 0, one finds from Eq. 
(I2.25> that the Ga{l) coupling decreases under scaling and is 
marginally irrelevant, while Gp{l) is marginally relevant. The 
relevance/irrelevance of the couplings reflects the low-energy 
properties having finite/zero excitation gap in the correspond- 
ing modes. This behavior correctly reflects the properties of 
the ID Mott insulator, where only the charge degrees of free- 
dom is frozen due to the finite Mott gap and the spin has gap- 
less excitations. By integrating out Eq. ( I2.25a> analytically 
from I — Oto I — Ip = ln(A/Ap), one can obtain the charac- 
teristic energy scale Ap as 

Ap = GpAVG;cxp(-l/2Gp), (2.26) 

where Cp is an integration constant depending on Gp{lp). 
This formula reproduces the exactly known Mott gap in one 
dimension in the weak U region, since the U dependence of 
Ap is given by Ap cx ytp7exp(-27rt|| /[/).'*** 

m. FORMULATION OF RENORMALIZATION GROUP 

In this section, we derive the RG equations for the QID 
half-filled Hubbard model in the two-loop level based on 
the Kadanoff- Wilson cutoff scaling scheme.' In the one-loop 




FIG. 3: The second-order diagrams contributing to the vertex cor- 
rections. The open square is the vertex for forward and backward 
scatterings, i.e., g\^, g2±, and and the shaded square is the 
one for umklapp scattering gs± and 33 y. The solid (dashed) Une 
refers to a right-moving (left-moving) electron, p — +( — ). The 
slashed line represents that the electron has energies in the shell 
Ai+di < |£p(fc)| < A; , while the crossed line represents the elec- 
tron having high energies determined by the momentum conserva- 
tion. The diagrams where the crossed line and slashed line are inter- 
changed are also taken into account. 

level, the formula tion for the QID case is found in Refs. 
'2r22'23'24'25'26'. In this scheme, we take partial integra- 
tion of the partition function over the fermion degrees of free- 
dom in the outer energy shell and scale the bandwidth cut- 
off A as A/ = Ae~' where I is the scaling parameter. We 
perform the logarithmic approximation, i.e., we keep the dia- 
grams which become logarithmic singular in the ID limit and 
thus the resultant QID RG equations are natural extensions 
to those for purely ID case. In order to simplify the nota- 
tions, we introduce the dimensionless couplings Gj^^q k^ j^^) = 
9v(qMM)l'^'^'"^ where v ^ p,a, c, s, cs. 

A. Peierls and Cooper bubbles in the one-loop level 

First we focus on the one-loop contributions due to the 
second-order vertex corrections. Possible Peierls and Cooper 
bubble contributions, due to the normal and umklapp scatter- 
ing, are shown in Fig.|3l We integrate out the electron degrees 
of freedom which have energy in the shell Ki+di < \£p{l^)\ < 
A;. The respective Peierls and Cooper bubbles have the trans- 
verse-momentum (i.e., patch-index) dependence of the exter- 
nal variables, as discussed in the literature. ^'•^^•^■'■^"' This effect 
is crucial to induce the transverse-momentum dependence of 
the coupling constants. There remain ambiguities in the se- 
lection of the longitudinal momenta for the external variable, 
since, in general, all the momenta of vertex cannot be set on 
the Fermi surface if the Fermi surface is warped. In this 
paper, we set three of four external momenta being on the 
Fermi surface and the longitudinal momentum conservation 
for each vertex (even for the internal momenta) is also con- 
sidered. The choice of the external longitudinal momenta, in 
addition to the transverse momenta, affects on the internal mo- 
menta and also on the RG equations. To keep the symmetries 
discussed in the preceding section, we also take into account 
the different choice of three of four longitudinal momenta on 
the Fermi surface. The explicit form of the Peierls bubble 
is given by -(T/V) Y^l""- J2n So+ (k, iuJn)Qo- (k - q, iuj,,) 
where Gop{k, it^n) = [i^^n ~ ^p{k)]~^ is the Green's function 
for the noninteracting case. By taking summation of the Mat- 
subara frequency and by performing the outer-shell integral 
over constant energy, this Peierls bubble contribution is given 
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—I{Tr-q+ki+k2,k' .ki.k2) where we have used the particle-hole 
symmetry. 



FIG. 4: The logarithmic-singular second-order diagrams for the 
Green's function, contributing the self energy. The notations are the 
same as in Fig.|3| Other types of second-order diagrams are not log- 
arithmic singular even in the ID limit and can be neglected. 



by {2t:vNi_) ^Y.kkq,kMM) 
I(qMMM) is given in the T 

A 



(q,k,ki,k2) 



EE 



dl, where the cutoff function 
limit by 

e(A + M5,fc,fc,(0) 



2A + Mq,fe,fe.(0 



(3.1) 



The quantity Ag kk'{l) being the functions of t±{l) is given 
by 



Aq^k,k'il) = 2t±{l)[cosk + cos{k - q)] 

- 2t±_ (l) [cos k' + cos(fc' - (?)] . 



(3.2) 



The second term in the rhs of Aq^k.k'{l) appears due to the 
longitudinal momentum conservation. In the conventional ap- 
proach, this term has been neglected,^'-^^-^^'^"* however is cru- 
cial to reproduce the known RG equations in the two-leg lad- 
der system (Sec. llVl l. We note //^ 



{q,k,k\ ,/c2) 



= 1 in the ID limit 



(t_L — > 0). The Cooper bubble contribution is also calculated 
in a similar way and can be expressed, after some algebra, as 



B. Two-loop self-energy corrections 

To go beyond the one-loop RG theory, we have to take into 
account two-loop self-energy corrections based on the second- 
order perturbation. The Fermi surface deformation can be 
taken into account by considering these corrections and has 
been discussed intensively by Dusuel and Doufot,^*" based on 
the zero-temperature formalism. Here we perform the finite- 
temperature formalism and take the T — > limit at the fi- 
nal stage of the calculation. The second-order self-energy di- 
agrams are shown in Fig. 0] In the second-order perturba- 
tion, there are two types of corrections to the single-particle 
Green's function Q: One is the corrections to the wave- 
function renormalization factor while the other contributes to 
the renormalization of the velocity and the interchain hopping. 
In the present RG scheme, the renormalization factor can have 
a transverse momentum dependence. So we assume that the 
Green's function takes a form 



iLUn — v{pk\\ — kp) + 2^"]" COS A; 



(3.3) 



where = 



,L I— 



(= Zfc^). The explicit calculation of outer 



shell integration of the diagrams in Fig.|4]yields 



dl 



GR^{k,iuJn) = ^Oi?(^'*^") ~ ^];;^E^S(q,fe,fc') [Jo{q-k,k') - Jl{q,k,k')GoRik,iuJn)] 



(3.4) 



q^k' 



for the right-moving electrons. The second-order coupling constants contributing the self-energy corrections are put into a form: 



/o2 /-^2 

'^T.(q,k,k') 



k.k' 



^2±(q,k,k') 



^'\\(q,k,k') + 2 ^3±(«,fe, 2^3±(7r-(/+fc+fe',fe,7r-fc') 



^2 

1^2 r< n ^r<2 

+ 9^3||(<j,fe,7r-A:') ~ (^3\\{q.k,Tr~k')<-'3\\{Tr-q+k+k' ,k,TT-k') + ^ (^311 (jr-q+fc+fe' ,fe,7r-fc') ' 



(3.5) 



We note that the umklapp scattering with the same spins G311 
also has finite contributions which are absent in the ID limit. 
The quantities Jo(g,/c.fe') ™d Ji[q,k,k') denote the cutoff func- 
tions due to the warped Fermi surface, which are also de- 
termined by the quantity ^ ^/(Z) [Eq. J3.2H . These cutoff 
functions Jq and Ji take different forms depending on the re- 
lation between Aq^k,k'{l) and A: For \ Aq^k,k'{l)\ < A, these 
are given by 



Joiq.k.k') = 2Aln 

Jl{q,k,k') 



4A + Aq^k,k'il) 



4A- Aq,k^k'il) 

16A2 



(3.6a) 
(3.6b) 



For \Aq,k±, (1)1 > A, 



k,k') 



2 A In 



AA+\Aq,k,k'il)\ 
2A +\Aq^k,k'{l)\ 

2A 



sgn{Aq^k,k'{l)), 
(3.7a) 
(3.7b) 



2A 



4A+\Aq,k.k'il)\ 2A +\Aq,k.k'{l)\ 



There remain subtleties in the integral region of outer shell^i 
here we adopt the simplest shell integral following Ref.lJ. The 
scaling deviation terms'''* have been neglected. 

The self-energy corrections proportional to Qf^j^ {k, itOn) in 
Eq. ( 13. 4t contribute to the wave-function renormalization fac- 
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FIG. 5: The third-order diagrams contributing to the vertex correc- 
tions, which have an order 0{G^dl) in the ID limit. The notations 
are the same as in Fig.|3| Other types of third-order diagrams have 
an order OiG^df) and can be neglected. 

tor z^^ . The explicit RG equation of the wave-function renor- 
maUzation factor is given by 

^Inzfc = -i^^G]^(q,k,k')Ji(<i.kM')- (3.8) 

-L q.k' 

The self-energy corrections proportional to Jq in Eq. ( 13 .4> 
contribute to the renormalization of the velocity and the Fermi 
surface deformation. To simplify discussions in the present 
analysis, we neglect the velocity renormalization, since this 
effect would only yield quantitative changes. The Fermi sur- 
face deformation can be extracted from these second-order 
corrections. Since the Fermi surface is given by Eq. M.ll . 
the Fermi surface deformation thus corresponds to the renor- 
malization of the interchain hopping. By noting that the self- 
energy contributions to the interchain hopping should have 
transverse momentum dependence cos k, the RG equation of 
the renormalization for the interchain hopping is given by 

~ Ghi.k,k') Mi,k,k') cask. (3.9) 

-L q.k.k' 

The renormalization to higher-order interchain hopping has 
been neglected. 

C. Two-loop RG equations 

In order to complete the two-loop RG theory, one has 
to take into account the next-to-leading logarithmic contri- 
butions to the vertex part. The two-loop vertex correc- 
tions can be calculated in a similar way to that for the self- 
energy connection. The third-order diagrams with the next- 
to-leading logarithmic contributions are shown in Fig. |5] and 
yield the renormalization of the vertex as Gi(^q^ki,ko) ^ 
Zi(qMM)Gi(qMM)^ where i = 1_L, 2_L, ||, 3_L, and 3||. 
Other types of diagrams are of the order 0{G^dP) which are 



already taken into account in the one-loop level. As is well- 
known in the ID case, the RG is formulated by deriving the 
scaling equations for the "renormalized" coupling constants 
Gi{l) = GiZi{l)z'^{iy-^-* where z{l) is the wave-function 
renormalization factor. In the present QID RG, by keeping in 
mind that the vertex has a transverse-momentum dependence, 
the renormalized coupling constants are defined as 

Gi(q^ki,k2)i^) = G^q,ki,k2) Zt{q:ki,k2)i^) 

(3.10) 

for the normal scatterings (i — 1_L,2_L, ||), and 

Gi(q^ki,k2)i^) = G^g^kiM) ^l{q:kiM)i^) 

(3.11) 

for the umklapp scatterings {i — 3 _L , 3 1 1 ) . The wave-function 
renormalization factor z^^ comes from the rescaling of the 
electron field operator Even in the two-loop vertex correc- 
tions, the cutoff function due to the warping of the Fermi sur- 
face appears, which is given by 

J2{q+k":kiM;k',k") — ^^Jliq+k" -ki,k' ,k") 

+ l^'^liq+k" -k2,k' ,k")- (3.12) 

The cutoff function /, Jq, Ji, and J 2 are not universal and 
would take different forms depending on the RG formulation. 
The function / [Eq. j3.1H is not continuous as a function of 
Aq^k.ki (0 which would be due to the sharp cutoff of the band- 
width. This unphysical discontinuity of / affects the results 
of the numerical integration of the RG equations. In order to 
avoid this unphysical effect, we replace J by a smooth func- 
tion which reproduce the limiting behavior of Eq. ( 13. Il l for 
small and large Ag^k.ki (0- 

From the straightforward calculation of the diagrams in 
Fig.|5] we obtain the two-loop RG equations for Gij^(g ^3), 

G'2_L(g,/ci,fe2)' '-'l!(g,fei,fc2)' G'3J-((3,fci,fc2)' ™d '-'3||(q,fci,fc2)- We 
note that, if we set iVj^ — 2 and if we neglect the umklapp 

scattering G3j_(q^kiM) and G^ni^q^kiM)' our RG equations re- 
produce the two-loop RG equations obtained by Fabrizic^^ in 
the two-leg ladder system at away from half filling. By us- 
ing Eq. ( I2.10L we rewrite the RG equation in terms of Gp, 
Gcr, Gc, Gs, and Gcs- For the system with the spin-rotational 
SU(2) symmetry, the coupling constants satisfy the relations 
given by Eq. ( 12.151 1. The full RG equations in this case is given 
in the AppendixlHI For the particle-hole symmetric Hubbard 
model, the coupling constants also satisfy the relations ( I2.17t 
and M.22\ . By using all these relations, the RG equations with 
the SU(2) X SU(2) symmetry are extremely simplified. The 
complete two-loop RG equations for the coupling constants 
Gp(qMM) ™d G„(^g^kuk2) are given by 
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^j^^v{q,ki,k2) 2N± ^ ^ \P'-'^iq;ki,k2;k')-l-{q.k' ,ki,k2) Pv(q;ki,k2;k')-^(7T — q+ki+k2,k',ki,k2)\ 

~ ^jyT G„{q,ki,k2) X] [^l(9':fel,fc') Jliq'M.k') + Gj:(q',k2,k') Jl{q'M,k') 

~ 4/V2" ^''('■''l^''^) X! [^S(g',-fci+9,fe') '^l(g',-fci+<?,fc') + j.,-, Jl(q, ^_k2+q,k') 

+ 4Jy2" 51 { ~ '^:^G^(^T^_q_q/+ki+k2,kiM)]^l'{q~ki+k',q-k2+k';k',k';q') 

'^^i^(,'^ — q—q'-\-ki-]-k2 - ki ,k2)^u{q—k\+k'.q—k2+k':k'.k'\q' ) J* <-^2(q+fc' ^k2\k' ^k' ~q') 
~^ ^ y I [^'^(?-?',fcl-9',fc2-g') ~ Qi'G'y(7r-q-q'+fci+fe2,fcl-i3',fc2-g')] T''(*:i-*:'!*:2-fe';fel,fc2;>3') 

~ 2^'^i-^-g-q' + ki+k2.ki-q' ,k2-q')^v{ki-k' ,k2-k'[ki,k2-,q')"^'^'2ik'\ki,k2\k' ,k' -q'), (3.13) 

where v = p^a and the sign function Q^, is 8p = +1 and = —1. The index of the scaling parameter I in the coupling 
constants G^(^q j^^ }^^) is suppressed. The coupling constants for the self-energy corrections, Eq. j3.5> . can be rewritten in terms 
of Gp and G^, as 

2 2 1 2 3 

GT.(q,k,k') ~ Gp(q,k,k') + 2'-^p{q.k,k')Gp(Tr-q+k+k' .k,k') + 3G^(q^fe^fc/) — 2^a{q.k,k')G a(Tr-q+k+k' ,kM') ■ (3.14) 

The quantities a^, /3y, 7^, and (5^ (i^ = p, a) are defined as follows. The quantities ai, represent the one-loop Peierls bubble 
contributions given by 

'^p{q;ki,k2\k') = 2Gp(q fej^ fc/) G p(^q]~i + Gp(q fej^ fc/) G p(^T^ _ qj^k2 + k' ,k' M) ^ p{iT-q+ki+k' ,ki ,k') Gp(q^k',k2) 

+ GG„(qk-^ki-f Grj(q^k',k2) ~ '^Grj(q^ki ,k') Gcr(TT-q+k2+k' ,k' ,k2) ^ '^Grj(Tr-q+ki+k' ,ki,k') G„(^q^k'M) 
+ Gp(7r-g+fei+A:',fci,fc') Gp(7r-g+fc2+/c',fe',fe2) + 3Gcr(7r-g+fei +*;' ,fci Gcr(7r-g+fc2+fc',fc',fc2) i 

<^(7(q;ki,k2;k') = '^Gp(q^ki,k') Ga(q,k' ,k2) + 2GCT(q,fei ,fc') Gp(q,fe',fc2) ~ 4GCT(q,fei^fe') ^cr{q,k' ,k2) 

~ Gp^q k-^ ki '^ Go.(7r-g+fe2+fe',fe',fc2) ^ Go.(7r-g+fei +fe' ,fei ,fc') G p^q ^i k^'^ + 2Go-(q,fei ,fc') Gcr(7r-g+fc2+fe',fe',fc2) 
+ Gp(7r_g+fcj+/j' fcj fc/) G^^q ki k^^i + Gcr^q k-^ k'-^ Gp(7r-g+fc2+fe',fe',fc2) + 2GCT(7r-g+fci ,fei ,fc') ^cr(q,k' M) 
~ Gp(Tr-q+ki+k' ,ki,k') Grj(Tr-q+k2+k' ,k' M) ~ Go.(7r-5+fci ,fci ,fc') Gp(7r_q+fc2+fe',fc',fc2) 
~ 2Go-(7r-q+fci+A:',A;i,fc') Gcr(7r-g+fc2 ,fc2) ■ 

The quantities represent the one-loop Cooper bubble contributions: 

Pp(q\ki,k2\k') = Gp(g_j,2 + fe/ j.^ j./)Gp(q_j,j_|_j,/ fe/ j.^) + 3Go-(^_j.2+fc',fci,fe')^'^(9-fei+fe'^fc''fc2)' 
Pa(q;ki,k2\k') = G p{^q^k^j^k' M,k')G a{q-ki+k' ,k' M) ''^ ^ (y(q-k2+k' ,ki,k')G p{^q-ki+k' ,k' M) 
+ '^Gc[q^k2+k' ,ki,k')Ga[q-ki+k' ,k' M)- 

Finally the quantities 7^ and 5^ represent the two-loop vertex contributions: 

^p{qi,q2;ki,k2;q') = ^p(gi,fci,fci-g') ^p('J2,fc2,fe2-g') ^ 3GCT(gi,fci,fci-g') GCT(g2,fc2,fc2-g')' 

I<y(qi,q2\ki,k2\q') = ^ p(qi,ki,ki-q') G p[q^^k2,k2-q') ^ ^(riqi .ki ,ki-q') Ga-(^q^ ^k2.k2-q') J 

^P(qi,q2-ki,k2:.q') = G p(^q-^ ,ki-q') Gp(7r-(j2 -g'+2/£2 ,fe2 ,fc2 -g') + Gp(7r-gi -g'+2fei ,fci ,fci -g' ) G p(^q^^k2,k2-q') 

3Gg-(gj^ fcj /jj^g') Gg- ( ^ _ — g' + 2 fe2 , fc2 , fc2 ^ ) ^ Gg- ( TT — gi — g' + 2 1 , fe 1 , fe 1 — ) Gg- ( g2 , ^2 , ^2 ^ g' ) ' 

^cr(qi,q2;ki,k2;q') — p{qi,ki,ki~-q') Gp(7r-g2 — g'+2fe2 ,fc2 ,fc2 -g') ~ ^p(ir — gi -g'+2/£i ,fci -g') ^ p{q2,k2.k2-q') 

G£r(g^ jfcj jfej — g') GCT(7r — g2— g'+2/£2.fe2,fc2— g') '"^O' (l" — gi ^ + 2fel , fci , fc i — g' ) Gcr(g2 ,fc2 ,fc2— g') ■ 



We have only kept the marginal scattering processes. In the purely ID case, it is known that the Gcs term has ir- 
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relevant canonical dimension. In the present case, some 
of the Gcs{q,kiM) couplings have a marginal canonical di- 
mension, however, the RG equation for the Gcs{q,kiM) ^'^^^ 
not appear explicitly since the correction due to this term al- 
ways appears in a form {G^s{q,kiM)-'^cs{7r-q+kt+k2MM))^ 
which shows the same /-dependence of {~Gcr(q,k i.k2) + 
Ga{n-q+ki+k2,kiM))^ seen from Eqs. il.lSal and ( I2.22cl i. 
For N± = 8, e.g., the number of independent coupling con- 
stants reduces to 300 instead of 1440 for without assuming the 
symmetries. If the transverse momentum dependences of the 
coupling constants are neglected, the ID RG equations [Eq. 
(I2.25H are reproduced. 

From the numerical integration of the RG equations, we 
can estimate characteristic energy scales. Here we focus on 
the renormalized interchain hopping and the charge/spin exci- 
tation gaps. The effective renormahzed interchain hopping is 
given by 

if = Aexp(-?_L), (3.15) 

where the quantity l± is determined from t±{lj_) = A. In the 
noninteracting limit, the interchain hopping scales as t±{l) = 
t±eJ, then l± = ln(A/t±) and the effective interchain hop- 
ping trivially reduces to the bare interchain hopping i]^ — t± . 
This quantity characterizes the dimensional-crossover energy 
scale, below which the system cannot be regarded as a one- 
dimensional system any more. In addition, the Fermi surface 
deformation can be determined by this quantity. By noting the 
relation i2.ll . the deformed Fermi surface is given by 

iCff 

fc|?(fc_L) = fcF + 2^cosfcx. (3.16) 

V 

It is known that the Fermi-surface deformation comes only 
from the renormalization in the high-energy regime, since the 
coupling constants which appear in the rhs of Eq. j3.9t are the 
irrelevant couplings. 

In the present RG scheme, the information of the charge 
gap Ap and the spin gap A^r can be extracted by focusing on 
the combination of the coupling constants; 

^'^+ = 1^ G^(^k-k' ,k,k) ; (3.17) 
^ k,k' 

where v = p,a. This interpretation can be justified by noting 
that the uniform charge/spin susceptibility is determined by 
these quantities'^ and by the field-theoretical approach for the 
two-leg ladder (N±_ — 2) case as will be shown in Sec. II VI 
A typical scaling flow is shown in Fig.|6l where we have set 
N±_ = 8. As a reference, the scaling flow for the ID case is 
also shown. The charge coupling Gp+ shows similar behavior 
to that in the ID case, while the spin coupling Ga+ becomes 
relevant and have a finite fixed point value G*_ = — 1. We 
note that the magnitude of several coupling constants becomes 
large and exceed the unity under the scaling procedure for I > 
l±. By focusing on this scaling behavior of Gu+{1), we can 
estimate the magnitude of the excitation gaps by 

A, = Aexp(-/,) (3.18) 



1 
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ti(l) /A / 
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I 

FIG. 6: (Color online) The scaling flows of the coupling constants 
Gp+ (l) and Gcr+ (l) and the interchain hopping ti_ (/) /A for A'^j_ = 8 
with fixed U/t\\ = 2 and = 0.05. The case for t± = is 

shown by the dotted lines. 



where the quantity l^, is determined from |G'y+(Z,y)| = c 
where c is a numerical constant. In the present numerical cal- 
culations, we will set c = 0.7 and A ~ 2vkp. As seen in Eq. 
(I2.26> . these ambiguity simply affects on the numerical factor 
and our choice reproduce well the exact results of the Ap in 
the ID case.^'^ The interchain-hopping dependence of Ap, A^, 
and is shown in Fig.0 The charge gap is suppressed due 
to the interchain hopping but is always finite even when the in- 
terchain hopping exceeds the magnitude of the charge gap. In 
the present bipartite QID half-filled Hubbard model, we find 
that the charge gap is always finite for U > 0. This is contrast 
to the results obtained from the chain-DMFTii22i22i^ where 
the metal-insulator (Mott) transition has been suggested for fi- 
nite interchain hopping at T = 0. This difference would arise 
from the difference in the treatment of the Fermi-surface nest- 
ing of the system. In the present model, the Fermi surface is 
always nested perfectly even for the finite interchain hopping 
where the nesting vector is (tt, tt). In our approach, we fully 
take into account this effect, however in the chain-DMFT, the 
warping of the Fermi surface is not taken into account. We 
expect that the Fermi-surface nesting would play crucial roles 
in the QID Mott transition, since the ID Mott insulator itself 
is realized even in the small U region due to the commensu- 
rability effect, which would be sensitive to the Fermi-surface 
nesting. By means of the present QID RG scheme, the effect 
of the nesting deviation will be reported elsewhereii^ 



IV. TWO-LEG LADDER MODEL 

To indicate the validity of the two-loop RG equations ob- 
tained in the preceding section, we apply it to the two-leg 
Hubbard ladder model with a bipartite lattice. This model 
has been investigated by the RG method combined with the 
analytical field-theoretical method^^^^i^^iSi^ and by the nu- 
merical DMRG method,'*^ '*'* and it has been clarified that the 
spin-gapped insulating state called the ZJ-Mott phase is real- 
ized. Lin, Balents and Fisher obtained the highly symmetric 



11 



0.3 



0.1 





''' A, eff , 










i4aSfz___i — 





0.1 0.2 0.3 



FIG. 7: (Color online) The charge gap Ap, the spin gap A^, and the 
characteristic energy scale t'f*, as a function of tx/ty for N ^ = 8 
and U lt\\ = 2. The dashed line represents the magnitude of the bare 
interchain hopping. 

SO(8) Gross-Neveu model as an effective theory in the low- 
energy limit by using the fixed-point behavior of the one-loop 
RG analysis.''^ They further discussed finite-energy spectrum 
based on this effective theory, however, it is not clear that 
this high symmetry still holds at finite-energy scale. Actu- 
ally, the RG method allows us to study the characteristic en- 
ergy scales in addition to the fixed point behavior, however, 
the naive one-loop RG is not sufficient to estimate the exci- 
tation gaps, since the RG method breaks down at the scale 
corresponding to the largest gap, as mentioned before. A 
promising method is to derive an effective theory by tracing 
out the gapped modes based on the field-theoretical treatment. 
However, in the present two-loop RG, the excitation gaps in 
the respective modes can be estimated without following the 
tracing-out procedure. This is not so trivial if the respective 
modes are not independent. In this section, in order to check 
the validity of the present method, we consider the two-leg 
ladder system, which is a minimal model of the spin-charge 
coupled systems, and confirm that this two-loop RG theory 
reproduces results obtained by the DMRG method and further 
analyze the excitation properties in detail by combining the 
field-theoretical bosonization and fermionization method. 

The model can be obtained from Eq. j2.H by simply setting 
N \_ = 2. The possible values of the transverse momentum are 
fc^ = and TT. From the symmetry requirements [Eqs. ( I2.15> . 
(I2.17> . and ( I2.22H . the number of independent coupling con- 
stants reduces to 8 instead of 30 for without assuming the sym- 
metries. To respect these symmetries and to make the physical 
picture transparent, we derive the effective low-energy theory 
by applying the bosonization and refermionization3i2iS 

First we apply the conventional Abelian bosonization to the 
Hamiltonian. The field operators of the right and left-moving 
electrons are written as 

'ipp,s,c{^) = — exp {ipkp.cx + ip fp^s.c) , (4.1) 
V27ra 



where p = +/— represents the right/left moving electron, 
s represents the spin, ( represents the band index: ( — 
+(-) for k± = 0(7r), and fc^.i = (7r/2 ± 2t^/v) [see 
Eq. ( 12. 7H . The technical details can be found in Refs. 

and .42, The chiral bosons obey the commutation re- 
lations [ipp^sxi^)^y^p-s'X'i^')] = ipT^ sgn{x - x') Ssy Si^x' 
and [(p+,s.Q,'P-,s'.c'] — ™^s,s' Sc^X'- The Klein factors ris^, 
which satisfy {?ys,Ci '7s',C'} ~ 2(5s^s"^c C'' ^'"^ introduced in or- 
der to retain the correct anticommutation relation of the field 
operators between the different spin and band index. To ex- 
press the electron fields in terms of the bosonic fields repre- 
senting physical modes, we define a new set of chiral bosonic 
fields (j)^^, (tfp_, (g^, and by 

<^p,,,C ^ <^p+ + C^p~ + + <€- , (4.2) 

where s =t/i= +/^- The commutation relations for these 
bosonic fields are [(j)P^{x),(j)^,j,,{x')] ~ ip{n/4)sgn{x — 

X')5^,y5r,r' and [(t)+^{x),(j)~,^,{x')] = i{n/4)S„^ty'Sr,r'- 

From Eq. J4.U the density operator is given by 

The convention of the Klein factors is the same as Ref. 
From this relation, one finds that the boson fields (/)p± can be 
interpreted to denote the "charge" degrees of freedom, while 
(j>a-± to denote the "spin" degrees of freedom. 

To appreciate two SU(2) symmetries in the effective theory, 
we next fermionize the (l)a+, 0<t-, and (f>p+ bosonic fields by 
introducing the Majorana fermions (n = 1, • • • ,6 and p = 
R/L = +/-): 

± ifp +^ep) ^ ^ exp {^p 20^+) , (4.4a) 

± + ^ep) ^ ^ exp {^p 20^) , (4.4b) 
\/ Z V Zna 

i= + ,^5) ^ ^ , (4.4c) 
\/ Z V Zna 

where is the Klein factor, satisfying {n^ri ^v'r'} — 
^v.v'^r.r' and = 1. These Majorana fields satisfy 

the anticommutation relations: {Cp (2^)7 Cp' (2^')} — ^ 
x') Sp^pi Sn.n' ■ The Hamiltonian can be refermionized in terms 
of the Majorana fermions. Our new finding is that the two sets 
of three Majorana fields form triplets, due to the constraint of 
two SU(2) symmetries. So we define 

— ("apiCpjCp); Cp = (CpjCpiCp)- (4.5) 

The g -olog y Hamiltonian / dxHuddav = (-f^o + Hi)\n^=2 
[Eqs. M.5\ and MM with N±_ = 2] can be reexpressed in a 
highly symmetric form as 
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Hi. 



adder 



-« 2 ' ^^^^ ^ ' 



9p+ 



5p- 



^gg-(0,0,7r) 

27ra 



27ra 

(li? ■ II) cos(20p_ + 8tix/w) 



27ra 



(Cfl •Cl)cos20p 



*5p(0,0,7r) 

27ra 



(4.6) 



where gp± = 5(5p(o,o,o) ±5p(77,o,o)) andg^i = i(.g^(o,o,o) ± 
.9cr(7r,o,o))- We note that the coupling constants and 
are the same as defined in Eq. ( I3.17t . From Eq. ( 14. 6> . one eas- 
ily finds that the 6 Majorana fermions are not independent and 
are grouped into two triplets ^ and In the derivation of the 
above effective theory, we do not use any fixed point values of 
the coupling constants but simply have used symmetry con- 
straints. This means that the structure of the theory maintains 
at finite energy scale. The physical meanings of the respective 
triplets becomes clear by noting the following relations. The 
total spin operator S [Eq. ( I2.13H can be expressed in terms of 
the Majorana fermions in a local form as S' = J dx J{x) with 



jy{x) 
r{x) 



(4.7) 



Similarly the total "charge" operator Q [Eq. 
expressed as Q = J dx J'{x) with 

j'y{x) = -^{eu.eB+^^eL). 



(4.8) 



up to the Klein factor Thus we find that the system has the 
"charge-triplet" excitations described by the C,p = (^p, Cp: Cp) 
Majorana fermions. The derivation of these relations is quite 
similar to that for the spin chains.''" These current operators 
satisfy the SU(2) Kac-Moody algebra at level k — 2.'*^ 

For the relevant interchain hopping, we also find high sym- 
metry in the p~ mode. In this case, the terms 5p(o,o.ir) 
and gcr(o.o,7r) in Eq. ( 14. 6> can be neglected due to the pres- 
ence of 8t±x/v in the cosine potential and then the effec- 
tive theory becomes SO(3)xSO(3)xU(l) symmetric, where 
the SO(3)xSO(3) is due to the formation of two Majorana 
triplets and the U(l) is due to the absence of the potential for 
the bosonic field (/)p_ . This picture is only valid for large in- 
terchain hopping, since the U(l) symmetry is retained dynam- 
ically while the SO(3) x SO(3) has a microscopic origin. 

The U and t± dependences of the charge and spin gaps and 
of the crossover energy scale are shown in Fig. |8] The 
U/t\\ dependence of the spin gap reproduce qualitatively the 
DMRG numerical results,'*^ while our RG approach would 
overestimate the magnitude of the spin gap. As easily seen 
from Fig. |8] the energy scales of the charge and spin exci- 
tation gaps are different in the whole region of U which 



is contrast to the analysis based on the one-loop fixed-point 
behavior^** 

Next we examine the fixed-point behavior of the present 
analysis. The fixed point values are 

3p(7r,0,0) — 9p{Tr,0,n) = ~.9cr(0,0,0) = 9<T{Tr,0,TT) = +9*, 
3p(0,0,0) = 5p(0,0,Tr) ~ .9cr(0,0,7r) = 3(T(7r,0,0) = 0, 

where we find g* / {2Trv) = 2 in the present case. This implies 
that the symmetry is dynamically extended in the low-energy 
limit. The effective theory in the low-energy limit has been 
analyzed in the one-loop RG scheme and is known to be de- 
scribed as the SO(8) Gross-Neveu model. This effective the- 
ory can easily be reproduced from Eq. (14. 6> . To this end, we 
fermionize the boson fields by introducing another set of 
Majorana fermions: 



1 



-el) ^ 



71 



exp(z20^_), (4.9a) 
exp(i20^_), (4.9b) 



where Kp_ is the Klein factor. These Majorana fields satisfy 
the same anticommutation relations as before. By using the 
Majorana fields ^" for n = 1, • • • ,8 and by inserting the 
fixed-point values into Eq. ( 14.61 1. the fixed point Hamiltonian 
can be expressed as 



'T/Cff 

rt 



laddc 




(4.10) 



which is called the SO(8) Gross-Neveu modeli^ Here we 
note that this symmetry enlargement occurs in the low-energy 
limit, where all the excitations can be regarded to have the 
same magnitude of the excitation gap. In the finite en- 
ergy scale, however, this symmetry does not hold and has 
SO(3)xSO(3)xU(l) as seen in Eq. for relevant inter- 
chain hopping. 

Finally we examine the magnitude of the excitation gaps 
for the remaining modes, and and we show how the 
low-energy effective theory in the small interchain hopping 
t ^ Ap can be described and how the trivial limit of ^ 
can be reproduced in this Majorana-fermion description. The 
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FIG. 8: (Color online) The excitation gaps, Ap (the charge gap). 
Act (the spin gap), Ay, Ag (the gaps in the Majorana fermion and 
see text), and the characteristic energy scale tf for N± = 2. 
(a) The (7/t|| dependence with fixed t±/t\\ = 0.1 and (b) the tx/ty 
dependence with fixed t//f|| = 2. The dashed line represents the 
magnitude of the bare interchain hopping. 

form of the Hamiltonian (I4.6> is vaUd even in the small t± 
region, however, the physical picture in the t± ^ limit is 
not so trivial. In terms of the Majorana fermions ^" {n = 
1 , • • • ,8), the Hamiltonian ( 14. 6> in the t± ^ limit can be 
rewritten as Hc«\t^^Q = 'H'^s + with 

n=4, 5,6,7 

+ f(eAei + e^ei + el.ei + ekI)', (4.11a) 

n=l,2,3,8 

- f (cki + ele! + eici - cki)' , (4.11b) 

where becomes relevant and g^r becomes iiTelevant. Here 
we adopt the notation CnSm which denotes n massless boson 
modes in the charge sector and m massless boson modes in the 
spin sector^" If one assigns that the bosonic phase variables 
'?^p± and (/)^_|_ describe the "charge" and "spin" modes respec- 



tively, the tj^ ^ limit may be interpreted as C^S^, where 
the gapless "spin" mode is described by the ^ = (^^, C^) 
fermion (the central charge is c = |) and the gapless "charge" 
mode is by the fermion (the central charge is c = i). The 
total central charge is consistent with that for two isolated 
Mott insulating chains c = 2, however, this picture is not cor- 
rect obviously. The correct understanding in the t — > limit 
is that the low-energy state is described by (705*2 where the 
Majorana fermions ^'^ and should be regarded to describe 
the charge and spin degrees of freedom respectively. From this 
interpretation, we can expect that the magnitude of the gap in 
the Majorana fermions and ^® shows nontrivial behavior as 
a function of tj^, since one {^J) is gapped and the other (i^*) 
is gapless in the t±^ 0, while these form the multiplet and 
are transformed into the U(l) bosonic field 6'p_ in the large 
interchain hopping. In order to estimate the t±^ dependence of 
the gap in the Majorana fermions and of Eq. ( 14. 6> from 
the numerical integration of the RG equations, we consider 
the following combination of the coupling 



37 = i^lgpiT.fl^T,) + 9p(ofl,T^)J{^t^a/v) 

+ 5<T(7r,o,7r) - 9a(nfi,-^)J{^i^a/v)], (4.12a) 

+ 5(T(7r,o,7r) + 5(7(0,0, 7r)^(8t±a/i;)], (4.12b) 

where J{8t±a/v) is a cutoff function satisfying J{x) « 1 for 
a; ^ 1 and J{x) ~ for a; ^ 1. For relevant interchain 
hopping, we have 57 = 58 = ^{9a{w,o,w) + 5p(7r,o,7r)), which 
would reflect the low-energy property of the 6'p_ boson mode, 
and for t we have gj gp and gg, — > g„ reproducing 
the single-chain limit. The excitation gaps for these Majorana 
fermion are also shown in Fig. |8j where we have estimated 
by A„ = Ae"'" with G„(;„) = 0.7 {n = 7, 8). The ground 
state of the present two-leg ladder system is known to be the 
D-Mott phase for arbitrary t± > 0, however, as seen from 
Fig.lSja), the crossover from the ID-like Mott insulating state 
(having large charge gap and small spin gap) to the insula- 
tor of the ladder, which has SO(3)xSO(3)xU(l) symmetry, 
takes place at t± « A^jt^^o where and the excitation prop- 
erties for the Majorana fermion and ^® undergo consider- 
able changes. By increasing [Fig.|SJb)], the effective 
interchain hopping is suppressed extremely and the multi- 
plet of the and splits into two isolated Majorana modes 
where the low energy excitations are described by the Majo- 
rana triplet ^ = (^^, ^■^) as a lowest-energy mode and by 
the Majorana singlet as a second-lowest-energy mode. This 
picture reproduces the low-energy properties of the Heisen- 
berg spin ladder systemsi^ 

The present estimations of the excitation gaps are also jus- 
tified by noting that it reproduces the known quantum critical 
behavior obtained in the extended Hubbard model including 
the intersite Coulomb repulsion. The detailed estimation of 
the extended Hubbard model is given in the Appendixicl 
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V. CONCLUSIONS 

In the present paper, we have derived the two-loop RG 
equations for the half-filled bipartite QID Hubbard model 
with the nonperturbative treatment of the interchain hopping, 
based on the conventional Kadanoff- Wilson approach. By 
considering finite number of ID chains we have treated the 
transverse momentum k± as the patch index and have ob- 
tained the RG equations which can be extremely simplified 
reflecting the symmetry requirements of the Hubbard model. 
By solving these RG equations numerically, we have esti- 
mated the magnitude of the charge and spin gaps and clari- 
fied that the charge gap is suppressed due to the interchain 
hopping but is always finite even when the interchain hopping 
exceeds the magnitude of the charge gap. In order to justify 
the present approach, we have analyzed the RG scaling flows 
in the two-leg Hubbard case {Nj_ = 2) in detail based on the 
field-theoretical Majorana-fermion description and have clari- 
fied that the low-energy excitations have SO(3)xSO(3)xU(l) 
symmetry for large interchain hopping. 
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By applying this transformation to the g-ology Hamiltonian 
(12. 8> . we find that the transformed Hamiltonian is given in 
the same form of Eq. (12. 8t . but the coupling constants are ex- 
changed as 

9l±{q,ki,k2) ^ ~.93±(g,fci,7r-fc2)j (A3a) 
92±{q,ki,k2) ^ -.92±(7r-?+fci+fc2,fel,fe2)i (A3b) 

while 5||(,,fei,fc2) and .g3||(q, fci,fc2)J ^''e unchanged. In the spin 
part there are constraints (I2.15> due to the spin-rotational 
SU(2) symmetry. By applying the duality relation (IA3> to Eq. 
( I2.15> . we can derive the pseudospin SU(2) constraints ( I2.22> . 



APPENDIX B: FULL RG EQUATIONS FOR THE 
SPIN-ROTATIONAL INVARL4NT CASE 



In this section, the full two-loop RG equations are given 
in the case for the spin-rotational SU(2) symmetric case [Eq. 
(I2.14H . without assuming the pseudospin SU(2) symmetry 
[Eq. ( I2.20H . These RG equations are valid for the extended 
Hubbard model including additional spin-rotational symmet- 
ric interactions, e.g., intersite Coulomb repulsions. 

The RG equation for the interchain hopping is given by 
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APPENDLX A: CHARGE-SPIN DUALITY RELATION 

In this section, we derive the pseudospin SU(2) relations 
(I2.22> from the spin SU(2) relations (I2.15> by using the 
"charge-spin duality" transformation. It is well known that 
the Hubbard Hamiltonian (12. 1> is transformed to itself with 
U ^ —U, under the particle-hole transformation for the spin 

51,52 



down only, • i.e 



(-1) 



(Al) 



Since the density operators are transformed as (rtj,;.| + 
^ under this transformation, the charge 

and spin density operators are interchanged. In the Fourier 
space with the linearized dispersion, Eq. (EIJ IS rewritten as 

Cp,|(A;) ^ Cpj(fc), Cp^i{k) ^ cl ^{{p^T,7^) - k). (A2) 



d 



1 



^T,n{q,k,k') Jo{q,k,k') COS k 



q.k,k' 
\ ^ ^2 



~ 4/V3" Gl,n(q,k,k') Jo{q.k,k')'^O^f'' (^1) 

where the second-order coupling constants contributing the 
self-energy corrections are put into forms: 



^2 

"-^Sn(g,fe,fc') 
^2 



Gp(g,fe,fc') + 
2 Gl(q.k,k') 



' 3 G'^{q,k,k') 



(B2a) 



2 J./) 



2 Gc(q,k,k'}Gc{TT-q+k+k' ,k,k') 



(B2b) 



The cutoff function Jo{q,k.k') is given by Eq. ( I3.6a> . In gen- 
eral, the cutoff function for the umklapp scattering contribu- 
tions Jpj^^ J. J,, J takes a different form from that for the nor- 
mal scattering ones Jo(q.fc,fc'), however, if the system has the 
particle-hole symmetry, these become identical Jo(q.fc,fc') — 
J' 

o(<?,fe,fe')' 

The RG equations for the coupling constants without the 
assumption of the spin-rotational SU(2) symmetry are given 
in symbolic form as 
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— C 



k> -L q',k' 

~ AN^ ^PilMM) ^3{q,ki,k2,q' ,k')i 



(B3a) 



q',k' 



in-^'^al{q,kuk2,k') + XTTJ X] ^<^2(9,fci,fe2,9',fe') + (^1 ^ ^2) 



^cr(q,ki,k2) ^^^j "SjqMM^q' ,k')i 

+ 8jv2" "c2(g,fei,fe2,9',fe') + ((5''^l''^2) ^ (-?, TT - /;;2, TT - fcl)^ 

k' 

~ 4]V^ ^c(q,fei,fe2) '^c3(g,/ci,fe2,g',fc')5 



(B3b) 



(B3c) 



where S^i and S,^2 represent the one-loop Peierls/Cooper bubble contributions and the two-loop (third-order) vertex contribu- 
tions, respectively, and S3 and Scs represent the two-loop (second-order) self-energy contributions. The respective terms are 
given explicitly in the following. The one-loop Peierls and Cooper bubble contributions are given by 



^pl{qMM,k') — -j\ppl^q^kx,k')Gp(q,k' M) '^G„(q^ki,k') Ga(q,k' ,k2)\hq,k' MM) 

~ 2 \pp(<l-k2+k' MM) G p{q-ki+k' M M) + 3^o-(g-fe2+fc',fei,fe') G'o-(g-fei+fe',fe',fe2)] -^C(g-fei-fe2,fc',fci 
+ 2[G 

c{qMM) Gc(qMM) + ^c(7r-g+fei+fe',fei,fe') G'c(7r+g-fe2-fe',7r-fe',7r-fe2) 
~ Gc(qMM) ^c(7r+g-fe2 — fe'jir-fe'.Tr-fe)] -^(-g,7r-fe',7r-fci,7r-fc2)' 

\Gp(^q^kiM) Gcr{qMM) ~ Gcr{qM,k') G „(^q^k' M)] I{q,k' M M) 

~ [Gp(q-k2+k' ,ki,k') + Gcr{q-k2+k' ,ki,k')] G fj(^q-ki+k' M M)-^C (q-ki-k2M MM) 
— 2[G'c(<j,fej,fc') Gc{q^k2M) ~ Gc{qMM) Gc{-lT+q-k2-k' ,TT-k' ,-!T-k2)]^{-q^Tr-k' ,7T-ki,7T-k2)' 
[Gp(q,kiM)Gc{qM,k2) - ^Gcr(q^ki,k') Gc(q.k'M) + "^G^^q^kiM) G c(T!~q+k' +k2M M)'\^(qM ,kiM) 
+ [Gp(^_g+fei+fe2,fei,fe') Gc(q-ki+k',k'M) + G„(Tr-q+ki+k2M,k') Gc(q-ki+k' ,k' M)]^{-K-q+ki+k2M 



^al{qMMM) 



^cl{q,k-i,k2,k') 



k2) 



(B4a) 



(B4b) 



fel,fe2) ■ 

(B4c) 



Due to the particle-hole symmetry of the present model on the bipartite lattice, the respective cutoff functions sat- 



isfy I'iqMMM 



I" 

{qM MM) 



^{qM MM)^ ^{—q^'^—k' 



-ki^7T—k2) 



'-{qM MM) 



, and /, 



C{q—k\—k2M MM) 
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I(Tr-q+ki+k2,k',kuk2)^ where I[q^k' MM) is given in Eq. The two-loop vertex contributions are given by 



i,k2,q' M) ^ p{q+q' ,k\,k2) p(q — ki+k' ,k' M — q') ^ p(q — k2+k' ,k' ,k' —q') 



+ 3Gcr(g-fei+fe',fe',fe'-q') Ga{q-k2+k' ,k' M -q')] '^2{q+k' M M;k' M -q') 
— 2G'p(Tr-g-g'+/£i+fc2,fci,fc2) [Gc{q-ki+k' ,k' M -q') G c(q-k2+k' M M -q') 
+ G'c(^_g_q'+fcj-|-fe'_fc'^A;'_q') G(,(TT-q~q' +k2+k' M ,k' -q') 

^ Gc{q-ki+k'.k'.k'-q') G c{-n-q-q' + k2+k' .k' M -q')\ •^2(q+k' MMM M -q') 
+ G p(^q—q\ki—q' ,k2—q') [Gp(fei — fe' ,fci ,fei — g') Gp^fe^ — fe' ,fe2 .^2 — g') 

+ 3G(j(i:j_i:'_fcj^_fcj^„g') Ga{k2-k' MM-q')\ ^2(-fe';-fei,-fc2;ir-fe',7r-*:'+g') 
2G'p(7r+g+g'— fei— fe2,Tr — fei+g',7r — fc2+g') [Gc(fci — fc' ,fei — g') ^'(.(fcj — fc' ,fc2 ,/c2 — g') 
~l~ G'c(7r — g' + fei+/c',fci,fci— g') G(;(7r — g'+/c2 + fe' ,^2 ,fe2 — g') 

^ Gc(ki-k' MM-q') G c{Tv-q' +k2 + k' MM-q')\ '^2(-fe';-fci,-fe2;7r-/c',7r-/c'+g')' 



(B5a) 



' — '(72{q,ki ,k2.q' M) 



Ga{q+q' MM) {p p(q-kx+k' M M ~q') G p{q-k2+k' M M ~q') 

~ G a{q-ki+k' .k' M -q') Ga(q-k2+k' M -k' -q')\ ■J2{q+k' MMM -k' -q') 
+ 2Gcr(7r-g-g'+/ci+fc2,fcl,fc2) G c{q-ki+k' M -k' -q' ) G c{TT-q-q' + k2 + k' M M -q') '^2(q+k' MMM M ~q') 
G a{q—q' M—q' ,k2—q') \G p(kx—k' ,ki,ki—q') Gpffej — fc' ,fc2 .^2 — g') 

^ Ga{ki-k' MM-q') Ga(k2-k' MM-q')\ '^2(-fe';-fei,-fc2;7r-fc',7r-fe'+g') 



+ 2Go-(7r+g+g' — fci — fc2,7r — fc2+g',7r — fei +g') G i,[ki—k\ki.ki—q') Gc(TT — q' + k2+k' ,k2.k2—q') 



^ '^2(-fc';-fei,-fc2;7r-fc',7r-fc'+g')' 



(B5b) 



"c2(g,fci ,fc2,g' 



,fc') — [2Gc(x_g_g' + fej+fc2^fcj^fc2) G, 



(7{q—ki-\-k' M M ~q') Ga(7T — q—q' +k2 + k\k' ,k' —q') 
G c{q+q' ,ki .k2) G p{^q—lzi-\-k' M M ^q') G p(^7r — q~q' +k2 + k' ,k' ,k' —q') 
+ Gc(q+q' MM) Ga{q-ki+k'M.k'-q') Go-(7r-g-g' +fe2+fe' ,fc' ,fc' -g')] J2(q+k' MMM M -q') 
[2Gcf7r — g — g'+A:i+/i:2 -fci — g' ,^2 — g') Go-(/i:i —k',k\.k\—q') Gfj(-^^q' — l^2 — k' .7T — k2.7T — k2+q') 
G c{q—q' ,ki—q' M—q') G p(ki —k' ,ki ,ki—q') Gp^jr+g' — fe2 — fe' — ^2 — fe2+g') 



+ G, 



c{q—q\ki~q'.k2—q') G —k',k\.k\~q') G fj(^^j^q' —k2—k'.'n — k2.TT — k2+q')^ 



y 7" 

^ '^2(-fc';-fci,-fe2;Tr-fe',7r-fe'+g')' 

and the two-loop self-energy contributions are given by 



(B5c) 



= G 



^2 T/ 

" ^Su(g',fci,fc') '^l(g',fci,/c') 
7' 

Y,u(q'MM) ■-^l(q',k2M) 



+ G|n(g',fc2,fc') ^Hq'MM) + 
+ ^ln(g',-fci+g,fc') '^l(g',-fci+g.fc') + ^Su(g',-fci+g,fc') Jl(q' -ki+qM) 
+ ^I;n(g',-fe2+g.fc') Jl(q' -k2+qM) + ^lu(g' ,-fc2+g,fc') "^((g' ,-fc2+g,fc') ' 
■^c3(g,fci,fc2,g',fc') = G|ii(g',fci,fc') + ^Su(g' ,fei ,fe') "^Ug' ,fci ,fc') 

+ ^I;n(g',7r-fc2,fc') Jl(q' ,-iT-k2,k') + ^Eu(g' ,7r-fc2 .fc' ) "^((g' ,7r-fc2 ,fc') 

+ ^ln(g',-fci+g,fc') '^l(g',-fci+g,fc') + ^lu(g',-fci+g,fc') '^l(g',-fci+g,fe') 



(B6a) 



+ g: 



En(g',7r+fc2-g,/c') Jl(q' ,ir+k2-qM) 



Su(g',7r+fc2-g,fc') '^l(g',7r+fc2-g,fc') 



(B6b) 



r 



The cutoff functions Ji and J( (J2, J2, and ) depend on 
the lattice geometry of the model and take different forms 
in general. However, in the present bipartite model, the 
respective cutoff functions satisfy t/^^ j, j,/) = Ji{q.kM)^ 
111 L 7 

'^2(qMM\k' M') ~ 'J2{qMM\k' M')- 



J' 

'^2{qMMM -k") 



We also obtain J 



1(— g,7r— fci ,-n — k2) 



J-. 



l(q,kiM) 



and 



'^2( — fc'; — fei , — fc2 ;7r — fc'.ir — fc'+g') 

particle-hole symmetric case. 



J2(k'MM\k'M-q') f^'" 



If the interaction is on-site one only, the system has the 
pseudospin SU(2) symmetry, where Eq. j2.22> is satisfied. 
By using Eq. ( I2.22> . the coupling constant G|,j^j.j,,^ = 

G|„(g,fe,fc') + fe^fc,) can be rewritten in terms of Gp and 
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Ga and reproduces Eq. J3.14> . Then the RG equation for 
the interchain hopping [Eq. (IBIH leads Eq. ( 13. 9t and those 
for the coupling constants [Eqs. ( IB3a> and ( IB3bH lead Eq. 
( I3.13> . The explicit RG equations for the umklapp scattering 
[Eq. (IB3cH can be suppressed due to the pseudospin SU(2) 
symmetry [Eq. J2.22H . 



APPENDIX C: EXTENDED TWO-LEG LADDER MODEL: 
CHECK OF QUANTUM CRITICAL BEHAVIOR 

In Sec.|in| we have estimated the magnitudes of charge and 
spin excitation gaps by using Eq. O.lSt . If the charge and 
spin modes of the system are decoupled, such as in the single 
chain case, this method trivially works since the coupling con- 
stants representing respective modes are decoupled. However 
in the present Ai"^-chain system where the charge and spin de- 
grees of freedom coupled with each other, one may consider 
that the present analysis does not work since the RG approach 
may break down at a energy scale corresponding to the largest 
excitation gap. In order to justify the present estimation of 
excitation gaps, we have considered the two-leg ladder model 
(A^j^ = 2) which is a minimal model with the spin and charge 
modes coupled. As akeady mentioned in Sec. II VI the U de- 
pendence of the spin gap [Fig. |8] (b)] shows similar behav- 
ior to the DMRG results'*^. In this section, we reconsider the 
two-leg ladder systems and we show another evidence which 
supports strongly the validity of the present estimation of ex- 
citation gaps. 

We consider a toy model including an additional interaction 
V' which denotes the next-nearest-neighbor Coulomb repul- 
sion. The spin mode in this model is known to exhibit quan- 
tum critical behavior within a nontrivial universality class. 
The purpose of the present section is to check whether the 
present method reproduces correct behavior of the quantum 
critical point (QCP). The Hamiltonian of this toy model is 
given by 

H' = (yc],l,s'^J + l,l,s + c\2,sCj+l,2,s + H.c.) 

j 

+ ^' X! ("j.i^-J+i^s + '^■J^2?^J+l,l) ■ (CI) 
J 

The notations are the same as in Eq. ( 12. H . This extended 
two-leg ladder model is examined by the field-theoretical 
method^°. For small V' , the rung-singlet (or D-Mott) state is 
realized where the ground state is unique. By increasing V' , 
this rung-singlet state changes into a spin-Peierls (or PDW) 
state (see Fig. 9 of Ref. Boh where the ground state has two- 
fold degeneracy and breaks translational invariance along the 
chain direction. From the field-theoretical approach, the quan- 
tum critical behavior is confirmed on the transition point be- 
tween the rung-singlet state and the spin-Peierls state. On this 




FIG. 9: (Color online) The V'/f y dependences of the charge gap Ap 
and the spin gap Ao-, for the extended two-leg ladder model {N±_ = 
2) with = 1 and fx/fy = 0.5. 



QCP, the spin gap collapses and the effective theory for low- 
energy states is known to be described by the c = 3/2 confor- 
mal field theory where c is the central charge. 

This extended Hubbard model can also be analyzed in the 
present framework of the two-loop RG, where the only dif- 
ferences from the analysis in Sec. [Ill] are that (i) the g-ology 
coupling constants in Eq. ( I2.8l l have explicit momentum de- 
pendence and (ii) the pseudospin SU(2) [Eq. ( I2.22H is not re- 
tained due to the presence of the additional interaction. The 
RG equations in this generalized case are given in the Ap- 
pendix [E] The estimated charge and spin gaps as a function 
of y/t|| is shown in Fig.|9] We find that the present approach 
reproduces the critical behavior since the spin gap becomes 
small around the QCP and collapses just on the QCP. The crit- 
ical value of V' is consistent with Fig. 9 in Ref. 40. The RG 
scaling flows on the QCP show that the coupling Gp+ reaches 
of the order unity for I > lp+, however, the coupling \ Gcr+ \ re- 
mains small and becomes irrelevant Gct+(oo) = 0. Such scal- 
ing behavior is the same as expected from the field-theoretical 
approach,"*" and thus the present results can be justified even 
for spin-charge coupled systems. 

From the technical point of view, we discuss the reason why 
the present analysis works even for spin-charge coupled sys- 
tems. If one of the coupling constants reaches of the order 
unity in the scaling flow, the RG method breaks down where 
it can be understood that the corresponding mode has an ex- 
citation gap. In order to analyze the lower-energy properties 
further, the gapped mode should be traced out and one should 
derive the effective low-energy theory for remaining modes. 
Then one can apply the RG method to it again. In this con- 
text, the quantum critical behavior was confirmed in Refs.*^ 
and 41. As for the two-leg ladder systems, we find that this 
tracing-out procedure almost corresponds to the replacement 
of the relevant coupling constants to unity. On the other hand, 
in the scaling flow of the present two-loop RG, the coupling 
constants remain finite even for the relevant ones (see Fig.|6jl. 
Thus one can consider that such trancing-out procedure of the 
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gapped mode is performed automatically in the present two- 
loop RG approach. The minor differences between these two 
approaches do not affect the numerical results. Thus we find 
that the present approach to estimate the different energy gaps 
works even for spin-charge coupled systems. 

Finally we note that the procedure of the derivation of the 



effective theory is not straightforward and restricted to the 
N± = 2 case only. In the present estimation based on the 
two-loop RG, there is no need to derive such low-energy ef- 
fective theory explicitly and thus this fact is the reason why 
it is easy to extend the analysis to the large number of chains 
systems. 
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